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ABSTRACT 

We have independently measured the genus topology of the temperature fluctuations in the 
cosmic microwave background seen by the Wilkinson Microwave Anisotropy Probe (WMAP). 
A genus analysis of the WMAP data indicates consistency with Gaussian random-phase initial 
conditions, as predicted by standard inflation. 

Subject headings: cosmology — cosmic microwave background: anisotropy 



1. Introduction 

The greatly anticipated results from the WMAP project (Bennett et al. 2003a) have redefined the state 
of the art in Cosmic Microwave Background (CMB) science. The WMAP project has provided full-sky cov- 
erage of the CMB at unprecedented angular resolution. The improvement over the only previous full-sky 
dataset (COBE, Smoot et al. 1992) is of order 500 in terms of resolution elements, and the sky cover- 
age improvement over balloon and ground-based experiments (e.g., de Bernardis et al. 2000 [Boomerang], 
Stompor et al. 2001 [Maxima], Mason et al. 2002 [CBI], Kovac et al. 2002 [DASI]) is several hundred. 
The experiment has resulted in dramatically improved constraints on cosmological parameters, such as the 
matter density, Q m , Hubble Constant, Ho, and cosmological constant, S7a (Spergel et al. 2003). 

The tool of choice for assessing the cosmological parameters, is the power spectrum, where one calcu- 
lates the products of the spherical harmonic coefficients a,£ m and their complex conjugates (for each I and m 
in the Yi m expansion of the CMB sky, summed over m values to form Cg). This ai m a\ m product, however, 
explicitly removes phase information. And while the power-spectrum is a very powerful tool for assessing 
cosmological parameters such as Q m and Q^, it does (explicitly) remove any phase information contained in 
the ai m coefficients themselves. These phases, however, contain critical information for characterizing the 
primordial density fluctuations. Namely, standard inflation (e.g., Guth 1981, Albrecht & Steinhardt 1982, 
Linde 1982, Linde 1983) predicts that the temperature fluctuations in the CMB, at the resolution measured 
by WMAP, will be characterized by spherical harmonic coefficients with Gaussian distributed amplitudes 
and random phases. The WMAP data provide our best opportunity to date to test that hypothesis. 

The genus topology method developed by Gott, Melott & Dickinson (1986) directly tests for the Gaus- 
sian random-phase nature of a density (or temperature) distribution in 3 dimensions (Adler 1981; Gott, 
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Melott & Dickinson 1986; Hamilton, Gott & Weinberg 1986; Gott, Weinberg & Melott 1987), or in 2 di- 
mensions (Adler 1981; Melott et al. 1989). Coles (1988) independently developed an equivalent statistic in 
2 dimensions. The 2 dimensional case has been studied for a variety of cosmological datasets: on redshift 
slices (Park et al. 1992; Colley 1997), on sky maps (Gott et al. 1992; Park, Gott, & Choi 2001; Hoyle, 
Vogeley & Gott 2002); and on the CMB, in particular (Smoot et al. 1992; Kogut 1993; Kogut et al. 1996; 
Colley, Gott & Park 1996; Park, C-G. et al. 2001). 

The WMAP team has recently measured the genus of the WMAP sky (Komatsu et al. 2003), and 
demonstrated that the WMAP results are consistent with the Gaussian random-phase hypothesis. To do this, 
they carried out a large number of simulations of the CMB, in which the spherical harmonic coefficients 
were drawn from a Gaussian random-phase distribution. They then used their known beam profiles, to 
simulate the results in each frequency, and applied the KpO (Bennett et al. 2003b) mask, just as one would 
with the real dataset. For each of these simulations, they computed the the genus (as defined by Melott et 
al. [1989]) and compared it to the to that of the real dataset. They found that the real dataset does not depart 
significantly from the genus of the simulated Gaussian random-phase datasets, very much in agreement with 
the results predicted for WMAP by Park et al. (1998). That group further explored Gaussianity by testing 
other Minkowski functionals (Minkowski 1903) and the bispectrum (related to the three-point correlation 
function), and found constraints on Gaussianity similar to those provided by the genus study. 

We seek, first, to confirm this result using our own methods, but also to compare the genus measured 
from the data directly to the theoretical genus curve for structures on a sphere. Rather than simulating 
Gaussian random-phase realizations of the CMB, we compare directly to the theoretical prediction for the 
genus in two dimensions. Furthermore, we provide some details about the nuances of carrying out the genus 
calculation on the HEALPix map projection and the stereographic map projection. 

2. WMAP Observations 

The Wilkinson Microwave Anisotropy Probe (WMAP) 3 project team has released its one-year dataset 
in several different skymaps measured at varying frequency and angular resolution. Namely, the project has 
measured the anisotropy at frequencies of 23, 33, 41, 61, and 94 GHz at angular resolutions of 0.82, 0.62, 
0.49, 0.33, and 0.21 degrees (FWHM beamwidths), respectively. 

The WMAP team has released the data in a unique map projection called the HEALPix projection 
(Gorski et al. 2000). We have chosen to re-present the WMAP results in Fig. 1 for a few reasons. First, 
we will be studying the genus on stereographic projections (sections 6 and 7), so it is useful to see the full 
sky map in that projection. Second, we have chosen a much different color scheme that is more relevant to 
genus studies than the original color scheme used by the WMAP team. The WMAP color scheme is pretty, 
informative and impressive, but is not symmetric with respect to hot and cold spots. In our color scheme, the 
mean temperature contour is white. Higher temperatures than the mean are represented as linearly redder 
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with temperature. Lower temperatures than the mean are represented as linearly bluer with temperature. 
Therefore, the amount of blue ink per pixel on the page is linearly proportional to the coldness of a cold 
spot, and the amount of red ink per pixel is linearly proportional to the the hotness of a hot spot. The scale 
runs from — 200^K (solid blue) to +200^K (solid red), the same range as in the temperature maps produced 
by the WMAP team (Bennett et al. 2003a). Since the two-dimensional genus of a Gaussian random-phase 
field is all about the symmetry between hot spots and cold spots, we find this color map highly instructive 
for our purposes. 



The properties of the genus are well-known in three dimensions (3D) (Gott et al. 1986, Hamilton et al. 
1986, Gott et al. 1987), but require some explanation in the two-dimensional (2D) case (Melott et al. 1989), 
particularly in the case of the sphere (Gott et al. 1990). 



For 2D topology on a plane the 2D genus of a microwave background map is defined as (Melott et al. 
1989): 



So for v > (/ < 0.5) there are more hot spots than cold spots, while for v < (/ > 0.5) there 
are more cold spots than hot spots. The genus is also equal to the integral of the curvature around the 
temperature contour divided by 2ir. If we were to drive a truck around an isolated hot spot we would have 
to turn a total angle of 2tt as we complete an entire circuit around the hot spot. Driving a truck around 
an isolated cold spot, we would turn a total angle of 2ir with the opposite sign, with a negative turn angle 
defined as one that is a turn to the left when the hot region is on your right. Thus, in a temperature field 
that is divided into pixels we may define a pixel as hot if it is above the contour threshold and cold if it is 
below the contour threshold. Then the contour line is actually composed of a series of line segments with 
turns occurring at vertices in the pixel map. There is a certain turn occurring at each vertex (which, divided 
by 27r, is equal to the contribution to the genus). Around each vertex there are four pixels, the contribution 
to the total turn (and genus) at that vertex depends on how many of the four surrounding pixels are hot and 



3. Genus Topology on a Sphere 



92D = Number of hot spots — Number of cold spots 



(1) 



For a Gaussian random field, 

g 2 D oc isexp(-v 2 /2), 
where v is a parameter that measures the area fraction in the hot spots: 



(2) 




(3) 
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This is carried out by an identical algorithm to that in the CONTOUR2D program as described in 
Melott et al. (1989). 

Now we wish to define rigorously the 2D genus on a spherical surface. Here we will follow our 
previous derivation (Gott et al. 1990). The 2D genus is defined to be equal to minus the 3D genus of the hot 
regions confined within a thin spherical shell, measured in 3D. Use lead paint to paint the hot regions on the 
surface of a balloon — burst the balloon and you will have solid curved lead shapes that will have a certain 
3D genus — take the minus of this number and that will be the 2D genus as we will define it. 

Recall that we have defined the 3D genus as the number of (doughnut) holes minus the number of 
isolated regions. Our 3D genus is equal to the integral of the Gaussian curvature over the contour surface 
divided by — Air. Thus, a sphere in 3D has a genus of —1, because it is one isolated region. The Gaussian 
curvature is +(l/r 2 ) while the area of the sphere is 47rr 2 , so the integral of the Gaussian curvature over the 
sphere is Air, independent of the size of the sphere, and the genus is —1. A cube has a genus of —1 also, 
because it has the same topology as a sphere. In this case, the curvature of the surface entirely consists of 
8 delta functions at the eight vertices where three squares meet at a point introducing a contribution to the 
Gaussian curvature integral of tt/2 at each vertex, equal to the conical deficit angle at the vertex. (If we 
sandpapered off the corners and edges of the cube, we would produce at each vertex an octant of a sphere 
containing a contribution to the Gaussian curvature integral of l/8th that of a sphere or tt/2. The edges 
would be sanded off to quarter cylinders which have no Gaussian curvature and the faces would have no 
Gaussian curvature, so all the curvature is located at the vertices). A doughnut has a genus of 0, because 
it has one hole and consists of one isolated region. A sphere with two handles has a genus of +1 on our 
definition because it consists of one isolated region and has two holes. 

So using these definitions we can define the genus of a microwave background map on the celestial 
sphere. Suppose for example we have one hot spot in the north polar region, and the rest of the sphere is 
cold. Then the genus would be +1, because the hot spot cap is one isolated region. Now imagine the hot 
region covers the northern hemisphere while the southern hemisphere is cold. The genus would still be +1, 
because a hemispherical bowl is one isolated region. Suppose the hot region covers all of the sphere except 
for a cold spot in the south polar region. The genus would still be +1, because this would look like a sugar 
bowl with out any handles, which is also one isolated region. The topology in each case is identical since 
one can be deformed into the other, so the genus should be the same in all three cases. If the hot region 
ran around the equatorial regions and there were cold spots at the north and south polar regions, the 2D 
genus would be because in 3D the equatorial band (confined in a thin spherical shell) is a doughnut in 3D 
which has a 3D genus of 0. If the sphere had N isolated hot spots (polka dots) on a cold background, the 
2D genus would be N. If the sphere had N isolated cold spots on a hot background the genus would be 
2 — N. (Remember that 1 isolated hot spot on a cold background [a hot spot centered on the north pole; 
the rest cold] is the same topologically as an isolated cold spot centered on the south pole [the rest hot], 
because the latitude of the contour line can be simply moved downward to transform one into the other.) 
Suppose we have a hot spot covering the northern hemisphere with the southern hemisphere cold. This is a 
spherical cap of area 2tt steradians. Its 2D genus is +1 because in 3D it is one isolated region. Considered 
as a spherical cap (part of a thin spherical shell) this isolated region must have an integral of the Gaussian 
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curvature on its surface equal to Air by the Gauss-Bonet theorem. There are several contributions to this: 
there is a contribution of 2ir from the outside (spherical) cap surface, and a contribution of 2-7T from the 
inside (spherical) cap surface (a spherical shell has an inside and an outside boundary both spherical), and 
a net contribution of from the delta function curvatures on the inside and outside edges at the equator 
marking the boundaries of the inner and outer spherical caps — the contribution of the outer edge is positive, 
while the contribution of the inner edge is saddle-shaped and therefore negative. The equator is a geodesic, 
so we could drive a truck around the equator without turning; hence, the integral of the turning around the 
contour boundary is zero. Since we are circling an isolated hot spot, why is the total turning not 2tt as it is 
on the plane? Because we are circling a curved region, and if we parallel transport (without turning) a vector 
around the boundary it will suffer a deflection equal to the integral of the Gaussian curvature in the interior 
of the boundary, which in this case is 2ir, the value we expect for a closed curve in the plane. The vector 
must return to its initial position after circling, so the total of the turning integral and the parallel transport 
deflection due to circling a curved region must add to give 2ir. 

Now consider a hot spot centered at the north pole and whose boundary extends to latitude 30°N. This 
is a spherical cap with area tt. It is an isolated region so the 3D genus is —1, and the sum of the integral of 
the Gaussian curvature over the spherical cap in 3D must be 47r. The integral on the outside of the spherical 
cap is 7r, and on the inside of the spherical cap is also tt, and the boundary surface connecting the inner 
and outer spherical cap shells is locally flat (part of a cone), so the contribution from the inside and outside 
edges of the boundary must be 2ir, which must be equal to twice the turning integral for a truck driving 
on the spherical surface (since the total Gaussian integral, Air, is twice the 2ir value for circling an isolated 
boundary [turning plus deviation from parallel transport]). If we were to parallel transport a vector around 
this boundary, we would find a deflection of ir because that is the area in the spherical cap; to return to its 
original orientation, a rotation of 2ir, the truck must turn an additional angle of tt. Indeed, if you would 
like to drive counterclockwise around the sphere at latitude 30°N, you will have to keep your steering wheel 
turned to the left, since this not a geodesic, and the total of that turning integral is tt. Now, locally it is the 
turning integral that is being calculated from the formula for a Gaussian random field. So let us define the 
effective genus as 

92D,eff = 92D ~ 2f, (4) 
where / is the fraction of the area of the sphere in the hot spots. In our example, where the one hot spot 
was the entire northern hemisphere, g2D = 1, / = 1/2 and g2D,eff = 0. The value of g2D,eff is equal to the 
turning integral of a truck driving around all the individual components of the boundary contour (around all 
the hot spots) divided by 2ir; in this case it is because the truck driving around the equator drives straight 
ahead without turning. For the case where the hot spot is centered at the north pole and extends to latitude 
30°N, / = 1/4, so the value of 92D,eff = 1/2, because the turning integral in this case is ir. Now for a 
Gaussian random field on the sphere 

92D,eff oc ^exp(-z/ 2 / 2 ), (5) 

because the Gaussian random field locally behaves on the sphere as on the plane to produce this contribution 
to the turning integral. Thus, in comparing the WMAP data to the random-phase formula we will use g2D,eff 
defined rigorously as described above. 
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4. Genus in the HEALPix Projection 

The WMAP data are plotted using an unusual map projection of the sphere, called the HEALPix pro- 
jection (Gorski et al. 2000). The sphere is effectively projected onto a rhombic dodecahedron. (This semi- 
regular polyhedron is the cell for the face centered cubic crystal lattice, the polyhedron representing the 
boundary of the set of points closer to a particular atomic nucleus than to any other.) It has 12 faces, four 
diamonds that meet at the top like the four sides of a pyramid, four diamonds that circle the equator, and four 
that meet at the bottom like the sides of an inverted pyramid. If we were to set this rhombic dodecahedron up 
aligned with the cardinal directions we could label the 12 faces as follows: TN (top north), TE (top east), TS 
(top south), TW (top west), NE (northeast, on the equator), SE (southeast, on the equator), SW (southwest, 
on the equator), NW (northwest, on the equator), BN (bottom north), BE (bottom east), BS (bottom south), 
BW (bottom west). There are 6 vertices where 4 diamonds meet at a point (i.e., TN, TE, TS, TW all meet 
at the top, and similarly TN, NE, NW, BN meet at a vertex). There are 8 vertices where 3 diamonds meet at 
a point (e.g., TN, TE, NE meet at a vertex). Each diamond may be subdivided into diamond shaped pixels. 
This is done by factors of 4, subdividing each diamond into 4 diamonds, and repeating by similarly subdi- 
viding each diamond. (In a similar way we may produce a checkerboard of 64 square pixels by dividing a 
square first into 4 quadrant squares, dividing each of these into 4 sub-squares, and each of these sub-squares 
into 4 checkerboard squares). In the WMAP case these are all diamonds (rhombuses) instead of squares, but 
the topology is the same. The subdivisions are made so that all the pixels are equal area, and diagonal rows 
of pixels touching point to point have centers that lie on circles of latitude. With N successive subdivisions 
of the diamonds into 4 sub-diamonds, we have a total of 12 x A N pixels of equal area covering the sphere. 
For example, WMAP uses N = 9, where there are 3,145,728 diamond shaped pixels covering the sphere, 
each with dimensions of approximately 0.11° x 0.11°, which is adequate resolution to show the WMAP data 
which has an angular resolution of approximately 0.3°. To study the topology properly, one needs pixels 
that are at least 2.5 times smaller than the smoothing length. 

Each of these diamond pixels on the sphere has angles which would be complicated to calculate, so 
a version of CONTOUR2D for the sphere would seem to be complex and require many calculations to 
calculate the turning angles. However, we may simplify this greatly by using the topological invariance of 
92D to projection (first from the sphere to the rhombic dodecahedron, then from the rhombic dodecahedron 
to a cube). This cube has six faces, which we will label T (top), N (north), E (east), S (south), W (west), 
B (bottom). Each diamond of the rhombic dodecahedron is oriented with its short diagonal along an edge 
of the cube, and each of the 12 rhombuses is mapped into two triangles that appear on adjacent faces of the 
cube. For example the rhombus TN is divided onto two triangles, one of which is mapped onto the T map on 
the cube and the other is mapped onto the N face of the cube. In fact, the letter designation of each rhombus 
tells us which of the two cubic faces it is mapped onto. For example, the SW rhombus is divided into two 
triangles which are mapped onto the S and W faces of the cube. The T face of the cube thus is a square 
which is divided into 4 right triangles whose hypotenuses form the sides of the square, and which meet in 
the center; these 4 triangles represent the upper halves of the four rhombuses TN, TE, TS, TW which meet 
at the top. 

In such a projection we paint the hot regions with lead paint on the surface of the cube, and the 3D 
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genus of these lead shapes is the same as if they were painted onto the sphere, because they have only been 
distorted in shape, not changed in topology. The total integrated Gaussian curvature over the cube, which 
totals Att is concentrated entirely in 8 delta functions at the 8 vertices, where there is an angle deficit of 7r/2 
at each, as 3 squares instead of four meet at each of these 8 vertices. Each of the 12 rhombuses of the original 
dodecahedron is mapped onto two right triangles meeting at their hypotenuses along an edge of the cube; 
these can be flattened out to make a square, with square pixels. So we have equivalently, 12 square maps of 
4 N square pixels each. Thus, we can consider each rhombus as a square map with square pixels and we can 
calculate giD for each using CONTOUR2D. Vertices that are within the rhombus are included. Vertices that 
are along an edge with the next rhombus are shared according to a prescription where the vertices along the 
top and right edges of the rhombus are assigned to it when laid out properly, and the bottom and left edges 
are assigned to adjacent rhombuses, since two rhombuses meet along each edge. That accounts for all pixel 
vertices except for those at the corners of the rhombuses. Each rhombus has 2 corners where 4 rhombuses 
meet (such as TN, TE, TS, TW meet at the top) and 2 corners where 3 rhombuses meet at a vertex (such 
as TN, TE, NE) which occurs at a vertex of the cube where there is a tt/2 angle deficit in our cubic map 
and where 3 pixels only meet at a point. The corners where 4 rhombuses meet at a point are like ordinary 
vertices in the plane where 4 square pixels meet at a point, so this is handled with CONTOUR2D. The value 
of the contribution to the genus from that vertex must then be divided and shared equally among each of the 
rhombuses that meet there. The prescription for vertices where 3 pixels meet at a point (at the corner of the 
cube in our projection) is as follows: 

hot 3 cold genus = the contour does not intersect the vertex 

1 hot 2 cold genus = 1/4 added to hot rhombus, added to others 

2 hot 1 cold genus = 

3 hot cold genus = 1/12 added to each of the three rhombuses. 

If we were to smooth the edges of the cube by sanding, the edges would turn into quarter cylinders 
(with zero Gaussian curvature) and the 8 vertices (where three squares meet at a point) turn into tiny octants 
of a sphere each with an integrated Gaussian curvature of tt/2. The faces of the cube are flat so all the 
curvature is located in the 8 vertices. The edges of the three rhombuses meet at 120° angles in the spherical 
octant. So in the case 1 hot and 2 cold there is a turn of 60° there in the contour plus a curvature of 30° 
(equal to 1/3 of the integrated curvature over the octant of the sphere) so the genus contribution (turn plus 
parallel transport deviation = 90°) should correspond to +1/4 (or tt/2 divided by 2tx). This is then added 
to the hot rhombus as this is the hot spot we are circling. In the case 2 hot, 1 cold, there is a —60° turn, 
plus two regions of 30° curvature contributing to the parallel transport deviation, making a total contribution 
of to the genus. In the case 3 hot, cold, the vertex lies entirely within the hot region but it adds tt/2 
to the curvature integral within that region and this must be shared equally with all three rhombuses to add 
+1/12 to the genus in each. The total genus for the sphere is calculated by adding the genus from all 12 
rhombuses. From g2D we can calculate <?2D,e// by subtracting 2f. We can calculate / as simply the fraction 
of the diamond shaped pixels that are hot, since these are all equal area. Thus, from the WMAP data we can 
calculate the 2D topology and compare it with that expected from a Gaussian random-phase distribution. 
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5. The Genus of the CMB, Measured by WMAP 

The most direct data product from the WMAP team is the "internal linear combination" (ILC) map of 
the CMB, which is given in the HEALPix format. This ILC map uses the optimum linear combination of 
the skymaps at the different frequencies to remove the Galaxy and some other foregrounds (Bennett et al. 
2003b). This is the primary map distributed as the best rendering of the CMB and has a resolution (beam 
width) of 0.3° x 0.3°. Though this map presents a few problems for directly computing the genus (which 
we will address shortly), it is worth checking the genus on this map directly. 

The HEALPix projection provides a natural division of the dataset (sky) into 12 independent regions. 
We therefore measured the genus in each of these regions and computed the total for the whole sky, which 
we plot in Fig. 2. In this figure v is computed in terms of the pixel-wise temperature mean and standard 
deviation measured over the whole sky (method 1). This method would be equivalent to the area fraction 
method described by equation (3) if the pixel histogram of temperature were strictly Gaussian. Fig. 3 is 
identical, except that v is computed by area fraction in each individual HEALPix rhombus (method 2); 
see equation (3). Using the area fraction method measures more directly the random-phase nature of the 
distribution, since it separates that from any departures from a Gaussian histogram of the temperature (which 
could be measured directly). In method 1, we use v defined in terms of the standard deviation in temperature 
over the whole sky; in method 2, we define v by area fraction in each rhombus, thus treating each of the 12 
rhombuses completely independently. The two methods provide a contrast of data treatment, and yet, as we 
shall see, give essentially identical results. 

The inner errorbars in Figs. 2 and 3 are the standard deviation of the mean in the genus for the whole 
sky, estimated from those the 12 independent regions of the sky (adding the values of <?2D,e// obtained from 
each). The standard one-sigma errorbars are shown as the inner errorbars. The outer errorbars show the 
95.4% (Gaussian two-sigma) confidence interval for a Student's-t variable with 12 degrees of freedom, 
which we will explain shortly. 

The solid curves in Figs. 2 and 3 give the Gaussian random-phase genus curve, according to equation 
(2), with the best fit amplitude applied. The amplitudes, A, [g2D,eff = Avex.-p(— v 2 /2)], in Figs. 2 and 3 
are A = 3432 and A = 3657, respectively. These are by far the highest in any genus measurements of 
cosmological structure to date. As expected, the first value of A is lower than the second, because the power 
at large scales (quadrupole and octopole) lifts some regions (rhombuses) to higher temperature, congealing 
some hot spots therein, but lowers other regions (rhombuses), congealing some cold spots therein. Method 
2 sets a median in each rhombus separately and thus discovers more structures. 

The Student's-t formulation is necessary since we are estimating the standard error with the distribu- 
tion itself. Namely, if the true genus for the whole sky at each point on the curve were \x and the true 
standard deviation in the mean expected for that number of structures were a (we expect the distribution to 
be Gaussian for Gaussian random-phase fields [Gott et al. 1990]), then our estimator for is the arithmetic 
average estimated from the 12 rhombuses (the values from each rhombus are multiplied by twelve, then they 
are averaged), x, and our estimator for a is s/y/n — 1, where s is the root-mean-square difference from the 
mean (as usual). We expect that ^ is given exactly by equation (2) and that our best fit for the amplitude has 
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negligible error compared to any one point on the curve. The expression (x — fi)/ (s/y/n — 1) is defined to 
be a Student's-t variate (Lupton 1993). 

On average, we would expect only one of the 21 points in each of the genus curves in Figs. 2 and 
3 to miss the curve outside of the 95.4% confidence errorbar. In fact, not even one misses. From the 
properties of the Student's-t distribution with 12 d.o.f., one would expect to find that 33.7% (~ 7) of the 
21 points miss the curve by more than the inner errorbars. In both cases, the agreement is again better than 
expected, with only three misses in each figure. This better-than-expected agreement can easily arise by 
chance from correlations between the genus measurements at neighboring values of v, as shown by Colley 
(1997) (these correlations mean that the 21 data points are not completely independent — neighboring genus 
values measure many of the same structures). Nonetheless, the genus of temperature fluctuations in the 
CMB, as observed by WMAP, is certainly consistent with the Gaussian random-phase hypothesis. 

6. Genus on a Stereographic Map Projection 

In the ILC HEALPix data release, there are a couple issues that must be addressed. The first issue is 
foreground emission. Bennett et al. (2003b) discuss in quite a bit of detail the foregrounds in the WMAP 
results. They provide a very useful sky mask that should be applied to avoid foreground contaminants, the 
principal one of which is the Galaxy. In a second study, we have therefore used their KpO (most conservative) 
sky mask, which also removes some 200 point sources. Furthermore, we have restricted ourselves to Galactic 
latitudes |6| > 18°. Since the mask must remove much of the area near the Galactic Equator, the HEALPix 
projection, with its 4 equatorial diamonds, ceases to be ideal. Instead, we re-project the map, using the 
stereographic projection for both the Northern and Southern Galactic Hemispheres (as in Fig. 1, in which 
the \b\ > 18° cut occurs at 72.7% of the total radius in each hemisphere). 

We can calculate giD from the north and south galactic caps by using stereographic map projections of 
the caps. Stereographic map projections of a hemisphere are conformal and preserve angles. A stereographic 
projection of the north galactic cap is a circle, in which azimuthal angle in the map is the galactic longitude 
in the cap. Radius in the map, r = 2tan[(90° — b)/2] (See Fig. 1). The outside boundary is a circle of 
radius r& = 2 tan [(90° — 18°)/2]. The contours on the sphere can be approximated as spherical polygons — 
geodesic segments joined at vertices where turning on the sphere occurs. The stereographic projection 
preserves angles, so the turning that occurs at vertices is mapped properly. But, geodesies (great circles) are 
mapped into arcs of circles rather than straight lines, which turn on the map but not on the sphere so the total 
genus g2D on the sphere is not automatically calculable from the total turning that comes from circling the 
density contour on the map. 

So instead of calculating g2D for the sphere from the stereographic map directly, which would be 
complicated, let us instead calculate the relation between <?2D,e// on a spherical cap, and <?2r>map measured by 
the CONTOUR2D program on the flat stereographic map of the cap. [Of course once we have <?2D,e// we can 
calculate the contribution to giD due to the structures in the north galactic cap by adding 2f • Area cap / (47r)]. 
But it is really g2D,eff we are interested in anyway because it is what will be compared to the theoretical 
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formula. Thus, we can skip a step by going directly after <?2L>e//- 

To begin, let us excise the spherical cap from the sphere by cutting it off and separating it from the 
rest of the sphere. This will add boundary terms. Each hot spot (think curved lead sheet) that intersects the 
boundary of the cap will be closed off by a curve along the boundary that will make it an isolated hot spot. 
Think about these hot spots in 3D when the cap is excised from the rest of the sphere: any hot spot that hits 
the boundary will become isolated, floating in space, because beyond the boundary there is nothing. This 
spherical isolated cap will have a value of <72D,excised cap = <?2L> excised map; the total genus of the isolated flat 
stereographic map sitting isolated in the plane — think of flat lead sheets lying in a plane; because one can 
be continuously deformed into the other. We can calculate the terms contributing to 5 f 2D,excised cap- 

52Aexcised cap = (2vr) _1 £ (turns on spherical cap) 
+/ • Area cap /(27r) 
+(27r)" 1 (turns at cap boundary) 
+/ 6 [l-Area cap /(27r)]. 

The first term, (2tt)^ 1 ^ (turns on spherical cap), is equal to <?2D,e//,cap for the interior of the spherical cap. 
It is the total turns made driving on the sphere around all the contours within the spherical cap, and excludes 
any boundary terms produced by excising the spherical cap. Now g2D,eff, cap ^ u exp(—u 2 /2), so it is 
what we want to calculate to compare with the theoretical genus formula. The term / • Area cap / (27r), where 
Area cap is measured in steradians, is the integrated Gaussian curvature within the hot spots in the cap divided 
by 2ir, which contributes to the genus by the 3D argument presented above. The term (2tt)~ 1 J2 (turns at cap boundary) 
is the sum of the turns that are made when you are driving along a contour and you encounter the cap bound- 
ary and then must turn a sharp corner there to continue driving along the cap boundary to complete circling 
that hot spot as an isolated region. In general, if along the boundary there are N hot segments, there will be 
2N such turns, as we enter and then exit each of these N hot segments along the boundary. Each of these 
2N turns will be some angle between and 2tt. The final term, ff,[l — Area cap /(27r)], represents the total 
turning on the sphere that occurs driving along the N hot segments of the boundary of the cap. If the entire 
boundary were hot then driving around it would have a total turn of 2ir — Area cap . (For example, if this were 
a hemispherical cap the area of the cap would be 2ir and the total turn circling the boundary would be zero, 
because the boundary is the equator and it is a geodesic so no turning is required to circle the boundary). 
As the area of the cap goes to zero, the curvature within the cap becomes unimportant to parallel transport 
and the total turn required to circle the boundary approaches 2n. In general, the angle deflection suffered by 
parallel transport around the closed boundary curve is equal to the integral of the Gaussian curvature inside, 
which in this case is equal to the area of the cap. When we have completed the circuit, we have returned 
to where we have started so we have suffered a total rotation of 2ir, so the turning we do on the sphere 
circling the boundary plus the area of the cap is 2ir. Thus, turning = 2ir — Area cap . The number fb is the 
fraction of the boundary that is hot, and so it is this fraction of the boundary that we drive around. Dividing 
this by 2tt, gives the contribution to the genus derived by driving along the hot segments of the boundary: 
fb[l — Area cap / (2ir)]. Using the Copernican Principle, one would not expect the boundary to be special, so 
on average the fraction of the cap boundary that was hot should be equal to the fraction of the sphere that is 
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in hot spots. On average then, we expect = /. Substituting we find: 

52D,exci S ed cap = 52D, e //,cap + ~ + ^- ^ (turns at cap boundary) + / ^1 - ^r^^J ■ 0) 

S2D,excised cap — 52Ae//,ca P + / + y~ ^2 (turns at cap boundary). (8) 

Now let us measure the genus on the excised stereographic map (<?2_D,extised map), recalling that the excised 
stereographic map is topologically equivalent to the excised cap. 

52D,excised map = tt~ 5^ ( turns on stereographic map) + — (turns at map boundary) + f b . (9) 

The first term, (2tt)~ 1 (turns on stereographic map) = g2D,map> is the quantity measured by CONTOUR2D 
on the flat map. CONTOUR2D examines only the Cartesian vertices within the stereographic map region 
and does not include any boundary effects. There is no / • Area term because the map is flat and there is 
zero curvature inside the hot spots. The next term, (2tt)~ 1 (turns at map boundary), adds the turns that 
are taken when the contour hits the outer circular boundary of the map. If we use a stereographic projection, 
which is conformal and preserves angles, each and every turn at the boundary in the map will be equal that 
encountered on the sphere, so: 

— ^ (turns at map boundary) = — ^ (turns at cap boundary). (10) 

The final term, /&, is the total turn taken on the N hot segments of the boundary of the map divided by 2ir. 
If we circled the entire boundary we would have to turn a total of 2tt, because the map is flat and Euclidean 
geometry applies. The N hot segments cover a fraction /& of the boundary (galactic longitude is mapped 
onto azimuth in the stereographic map, so the fraction of the map boundary occupied by the N hot segments 
is exactly the same as on the sphere). Thus, the total turn is 2tt f^, and dividing this by 2ir gives /&. As we 
have remarked above, if the boundary is not special we expect /& = /. Since the map is excised and sits 
alone on the plane, these outer boundary segments must be included to create isolated hot spots in the plane 
out of the hot regions that hit the boundary. Equating <?2D,excised cap and <?2D,excised map (because the cap can 
be deformed into the map), and substituting from the above equations we find: 

52D,excised cap = 52D,excised map • (11) 

52Ae//,ca P + / + ^2 ( turns at ca P boundar y) = #2D,ma P + ^2 (turns at map boundary) + /. (12) 

<72D,e//,cap = ff2D,map- (13) 

So "what you see is what you get." If we measure the 2D topology in the interior of the stereographic map 
using CONTOUR2D (which ignores the boundary) our result will equal, on average, 52D,e//,cap which is 
what we were looking for. This works for any radius spherical cap, since terms proportional to the area of 
the cap cancel out. And since the map is conformal, the turns at the boundary are equal in both cases and 
cancel out as well. Finally, since we expect the boundary to be typical, its coverage with hot spots should 
be equal to the fraction / of the sphere covered by hot spots, and this causes other terms to cancel as well, 
leaving the above simple result. 
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7. Genus Topology on Stereographic WMAP Results 

We may now compute the 2-D genus on a stereographic re-projection of the WMAP data, in which we 
address the issues of foreground contamination and noise variability arising from the scanning strategy of 
the WMAP spacecraft. 

First, to alleviate contamination from the Galaxy, we completely exclude the equator region up to 
|6| < 18°. The KpO mask does include removal of local regions for point source avoidance, and some 
additional Galaxy masking above the \b\ > 18° cut, particularly near the Galactic Center, the Gum Nebula 
and Orion Complex. Any vertex in the stereographic map projection that falls within the mask is simply 
ignored by CONTOUR2D. 

We proceed by dividing the stereographic projection into octants which correspond to octants on the 
sphere. These octants form 8 independent regions in which to measure the genus. A slight adjustment for the 
actual solid angle in each octant is made, because the KpO mask changes the solid angle of each octant very 
slightly (a few percent variation). We normalize the genus to the mean solid angle of the masked octants, 
and proceed treating each masked octant as one of 8 independent measurements of the genus per average 
masked octant. 

In carrying out the re-projection to the stereographic map, we have applied a slight smoothing to a 
total smoothing FWHM of 0.35°. Figs. 4 and 5 show the genus for the masked stereographic projection, 
computed by methods 1 and 2 (same methods as in the HEALPix genus computation). The errorbars are 
computed now for 8 octants instead of 12 diamonds, which changes somewhat the relative size of the 95.4% 
errorbars and the direct one-sigma errorbars. The best fit theoretical genus curve (equation 2) is plotted in 
each figure, with respective amplitudes of A = 2057 and A = 2186 (as with the HEALPix genus curves, the 
first value is slightly less than the second, as expected). The amplitudes have come down somewhat from the 
HEALPix maps due to the exclusion of the equatorial band (so the total genus covers only a fraction of the 
whole sphere), the additional masked pixels in the KpO prescription, and the additional smoothing applied 
in the stereographic projection. 

In each of the figures, the agreement with the theoretical genus curve is striking. Quantitatively, we 
would expect on average that one of the 21 points sampling the curve would miss the theoretical fit outside 
is 95.4% errorbar. With method 1 (Fig. 4), no points miss the theoretical curve outside their 95.4% errorbar, 
while with method 2 (Fig. 5), two points miss by a small margin. The one-sigma (inner errorbars) for the 8 
degrees of freedom of the octants reside at the 65.3% confidence interval, so on average, one would expect 
approximately 7 of the 21 points on each of Fig. 4 and Fig. 5 to miss the theoretical curve outside their 
one-sigma errorbars. Using method 1, we find that 3 points miss at this margin, and 9 miss using method 
2. Though the method 2 results show slightly more than the average number of misses expected at both 
the one-sigma and 95.4% level, the excess is certainly within the uncertainty arising from small number 
statistics in each case. 

The second contaminant in the ILC HEALPix map is that the noise behavior is very complicated, and 
differs substantially pixel to pixel. First, the beams at each frequency in the combination have different 



-13- 



widths. Second, the scan method of the WMAP satellite yields many more samples at the ecliptic poles 
than at the ecliptic equator (see Bennett et al. 2003a, Fig. 3). To alleviate some of these problems, we 
have downloaded from the WMAP website maps in the same 5 frequencies as used in the ILC map, but 
ones that have been pre-smoothed to the same effective beam-width (0.82°). These are then added with the 
same weights as used in producing the ILC map. At cost of number of structures, these maps significantly 
enhance the signal-to-noise of remaining structures, thus reducing the significance of the noise differences 
arising from the scanning strategy. This is, therefore our most conservative map. 

In projecting these maps into stereographic maps, we have smoothed to 0.9°, slightly more than the 
0.82° smoothing of the pre-smoothed maps. Figs. 6 and 7 show the genus measured in these stereographic 
projections. The best-fit amplitudes have dropped to A = 772 for method 1 and A = 831 for method 2, due 
to the smoothing. 

From these final stereographic projections, the genus measured at the 21 points should also be expected 
to miss the theoretical curve outside of the 95.4% errorbars once on each figure. What we find is that with 
method 1 (Fig. 6), no points miss the curve outside the 95.4% errorbar, while with method 2, (Fig. 7), one 
point misses by a small amount. The points miss the curve outside their one-sigma errorbars 2 times for 
method 1, and 4 times for method 2, again a bit better than expected. The final figure, Fig. 7, represents 
our most conservative estimate of the true genus available in the WMAP distribution, and shows excellent 
agreement with the Gaussian random-phase hypothesis. 

8. Discussion 

As we have shown, topologically, within the errors, the CMB measured by WMAP appears to be 
consistent with Gaussian random-phase. Indeed, the overall visual impression is strikingly random-phase. 
The stereographic projection is conformal, so shapes are preserved locally, though the scale grows as one 
approaches the the outer boundary (by a factor of two from the pole to the equator, as seen in Fig. 1). We 
may ask if there are any features visible in the map at all that look non-random-phase in any way. 

Our choice of map projection has a miscellaneous advantage over the Molleweide projection used by 
the WMAP team, which is that it shows the Galactic polar regions very well. In the south polar region, 
there is one possibly non-random feature, a fairly narrow red (hot) feature stretching from approximately 
the 8:30 position toward the 2:30 position curving like a bow below the South Galactic Pole. This feature 
is more visually apparent when viewed in landscape orientation. It has the appearance of a scanning error 
approximately along a great circle, but is not coincident with any particular locus in the WMAP scanning 
geometry. 

We have considered several possibilities to explain such a feature. First, there is the simple phenomenon 
of the "canals" on Mars, the linear features reported in naked-eye observations of the red planet in the late 
19 th century. We now know that there are no such canals, and that our eyes' highly evolved ability to 
detect linear features or edges was simply creating an illusion. Most likely the same thing is going on here, 
particularly since the feature is more apparent when the figure is viewed in landscape rather than in portrait 
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orientation. In fact, if one had a true Gaussian random-phase distribution and there were no features that 
looked unusual to the eye, that would be unusual. So this feature may be merely coincidental. 

What are the other possibilities? It is not coincident with the Magellanic Stream measured seen in 
21cm (Brooks 2000). Furthermore, we would expect the stream to have similar microwave "color" to the 
gas in the Galaxy, and hence be subtracted about as well as the Galaxy is in the ILC map. To test this, we 
constructed a three color map, using the K band as red, Q as green and W as blue, and found that while the 
foregrounds were distinctly pink to orange to yellow, the CMB was quite grey, and so was the feature. 

Another possibility that comes to mind is cosmic strings. Cosmic strings, however, do not amplify 
temperature along a line. Instead, a moving cosmic string produces a cliff, hotter on one side (the trailing 
side) than the other side (the leading side). Thus a hot stripe would require two parallel cosmic strings, 
having just passed each other (an interesting theoretical possibility [Gott 1991], but not a likely one, since 
we are unlikely to find the two most prominent cosmic strings within the horizon in such an alignment by 
chance). Also, there is no obvious evidence for cosmic string cliffs anywhere else in the WMAP data, which 
we have tested for by making gradient maps as suggested by Gott et al. (1990). 

The feature could also be due to the integrated Sachs-Wolfe effect from a great wall seen edge-on, but 
in looking at several galaxy catalogs and X-ray maps, we could detect no obvious feature with appropriate 
geometry. 

This list of possibilities is by no means exhaustive, and, in the end, it may be most likely that we are 
simply falling prey to "Martian canals" in an overall Gaussian random-phase map. Since there are countless 
possibilities for what one may see and deem non-random (circles, squares, triangles, even letters or words), 
the a posteriori statistical significance of any such weak feature is difficult to estimate. 

The overall consistency with of the genus topology in the WMAP data with that predicted a priori 
from a Gaussian random-phase distribution is quite spectacular. WMAP's unprecedented combination of 
angular resolution and all-sky coverage of the Cosmic Microwave Background has presented by far the best 
confirmation to date of the standard inflationary prediction of Gaussian random-phase initial conditions. 

In standard inflation, the structures seen in the CMB are fossil remains of random quantum fluctuations 
occurring just 10~ 33 seconds after the Big Bang. In galaxy clustering surveys, these fossils have become 
somewhat distorted due to non-linear effects and biased galaxy formation. In the CMB, however, we are 
seeing these fluctuations while they are still in the linear regime, making the WMAP measurements of the 
CMB a very attractive dataset for studying these primordial perturbations. 

With the hundreds of structures observed even at the 0.9° scale, the WMAP results have provided 
the largest and best dataset for studying the primordial structures. The genus topology test, in particular, 
tests the dataset for the standard inflationary prediction that the fluctuations should derive from a Gaussian 
random-phase distribution. All of our tests have confirmed this prediction dramatically. 
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Fig. 1. — The WMAP data (internal linear combination) in the north (top) and south (bottom) Galactic 
hemispheres, in stereographic projection, with the Galactic Center at center. The color map is a linear scale, 
chosen such that hot spots are red and cold spots are blue with a total range in temperature of — 200/uK to 
+200 / uK. The mean temperature contour is in white. 
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Fig. 2. — Total genus, <72D,e//> f° r the original HEALPix projection of the WMAP data (internal linear 
combination) for the entire sky, with no masking. The abscissa, v, is the departure from the pixel-wise 
temperature mean for the whole sky in terms of pixel-wise temperature standard deviation for the whole 
sky. Inner errorbars are computed from the standard deviation of the 12 HEALPix diamonds; the outer 
errorbars reflect the sigma value scaled for 95.4% confidence in the Student's-t distribution. 
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Fig. 3. — As in Fig. 2, except the v values are computed by area fraction (see text). 
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Fig. 4. — Total genus, g2D,eff> f° r tne WMAP data (internal linear combination), projected stereographically 
to a total smoothing of 0.35°, then masked according to the KpO standard. The abscissa, v, is the departure 
from the pixel-wise temperature mean for the whole sky in terms of pixel-wise temperature standard devia- 
tion from the whole sky. Inner errorbars are computed from the standard deviation of the 8 octants; the outer 
errorbars reflect the sigma value scaled for 95.4% confidence in the Student's-t distribution. 
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Fig. 5. — As in Fig. 4, except the v values are computed by area fraction (see text). 
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Fig. 6. — Total genus, g2D,eff, for the pre-smoothed WMAP data, combined according to the internal linear 
combination coefficients, then projected stereographically to a total smoothing of 0.9°, then masked accord- 
ing to the KpO standard. The abscissa, v, is the departure from the pixel-wise temperature mean for the 
whole sky in terms of pixel-wise temperature standard deviation from the whole sky. Inner errorbars are 
computed from the standard deviation of the 8 octants; the outer errorbars reflect the sigma value scaled for 
95.4% confidence in the Student's-i distribution. 
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As in Fig. 6, except the v values are computed by area fraction (see text). 
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